Heat Transistor: Demonstration of Gate-Controlled Electron Refrigeration 
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We present experiments on a superconductor-normal metal electron refrigerator in a regime where 
single-electron charging effects are significant. The system functions as a heat transistor, i.e., the 
heat flux out from the normal metal island can be controlled with a gate voltage. A theoretical model 
developed within the framework of single-electron tunneling provides a full quantitative agreement 
with the experiment. This work serves as the first experimental observation of Coulombic control 
of heat transfer and, in particular, of refrigeration in a mesoscopic system. 
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In conventional transistors such as a field-effect or 
single-electron transistor [l| , the electric current or volt- 
age is controlled by a gate. More recently, a different 
kind of working principle was demonstrated in meso- 
scopic structures: the electron transport can also be con- 
trolled by manipulating the quasiparticle temperature, 
i.e., the energy distribution of charge carriers [3, 
In these experiments, such manipulation resulted typi- 
cally in an increase of electronic temperature. In subse- 
quent experiments, the supercurrent of a superconductor- 
normal metal-superconductor (SNS) Josephson junction 
was controlled not only by heating but also by cooling 
the electrons in the active region [J]. A typical electron 
refrigerator consists of a SNS system with two insulating 
barriers defining a SINIS tunnel structure [l| . Biasing the 
device with a voltage of approximately twice the super- 
conducting energy gap results in extraction of heat from 
the N island, and hence a drop in its temperature. Ultra- 
small structures operated at low temperatures display 
pronounced charging effects which may profoundly affect 
the heat transport dynamics in mesoscopic systems. 

In this Letter, we investigate a SINIS structure with 
very small tunnel junctions which behaves as a gate- 
controlled single-electron refrigerator. In such a system, 
we demonstrate experimentally gate modulation of the 
heat flux by more than a factor of three. Our experimen- 
tal findings are successfully explained through a model 
which associates the thermal current with single-electron 
tunneling [6( . This structure provides an ideal framework 
for the investigation of the interplay between charging ef- 
fects and heat transport at mesoscopic scale. 

Figure []Ja) shows the measured heat transistor with a 
sketch of the measurement setup. The core of the sample 
consists of an N island, tunnel- coupled to four supercon- 
ducting leads, and a gate electrode with capacitive cou- 
pling to the island. Like in an ordinary single-electron 
transistor, the charging energy penalty for tunneling elec- 
trons can be modulated by the gate voltage. As we shall 
show in the following, in the heat transistor under suit- 
able bias conditions, Coulomb blockade allows the mod- 
ulation of both the electric current and the heat flux. For 
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FIG. 1: (color online) (a) Scanning electron micrograph of the 
heat transistor, showing also the measurement setup. In the 
middle, the normal metal (Cu) island subjected to cooling. 
Four superconducting probes (Al) with Al-Ox tunnel junc- 
tions can be seen above, and truncated capacitor plates (Al) 
with direct NS contacts on the left and right, (b) Theoreti- 
cal heat flow from the N island at the optimal bias point at 
different temperatures as a function of the gate voltage. T 
denotes the common temperature of the island electrons and 
quasiparticles in the superconductor, and T c = A/(1.764A;b) 
is the superconducting critical temperature. 



a direct verification of the Coulomb blockade effect, we 
also fabricated and measured a reference sample with the 
same geometry and parameters, except for large capaci- 
tor pads (area 5 x 10 4 /xm 2 ) attached to the island with 
NS contacts, thus suppressing the charging energy but 
ideally maintaining the thermal isolation of the N island. 
In the actual heat transistor we preserved only small su- 
perconducting protrusions instead of the large pads in 
order to maintain the same topology with a minimal ef- 
fect on the charging energy. The samples were fabricated 
with electron beam lithography and three-angle evapo- 
ration on an oxidized silicon wafer. The normal metal 
island consists of a copper block with lateral dimensions 
180 nm x 2300 nm and thickness 20 ± 3 nm. The super- 
conductor metal was aluminium, and the Al-Ox tunnel 
junctions had an approximate area 120 x 40 nm 2 . In the 
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heat transistor, the two inner junctions with 110 k£! re- 
sistances were used as a voltage-biased cooler pair. Con- 
sistent with the transistor terminology, we denote the 
bias voltage over the cooler junctions Vbs- The outer 
junctions with resistances of approximately 150 kf2 each 
were current-biased at picoampere level, enabling accu- 
rate thermometry of the island electrons by exploiting 
the temperature sensitivity of a NIS junction I-V curve 
in the 100-500 mK range. As only the DC properties of 
the device were relevant, all signal lines to the sample 
stage were low-pass filtered. 

Figure QJb) illustrates the expected performance of 
the transistor, displaying the modulation of heat flow 
(Q) with gate voltage (V g ) at different operating tem- 
peratures. The electric and thermal transport properties 
of a Coulomb blockaded SINIS structure can be derived 
through the theory of single-electron tunneling assuming 
no exchange of energy due to environmental fluctuations. 
To obtain a complete model of the experimental setup, we 
consider a normal metal island with four tunnel junctions 
to superconducting leads numbered 1 ... 4 with junction 
resistances and capacitances given by Ri and Ci, respec- 
tively, and the electric potential of each lead denoted by 
Vi. The gate electrode is coupled to the island with ca- 
pacitance C g . The total capacitance of the normal metal 
island reads Cs = C g + J2t=o ^ > wnere Co denotes the 
sclf-capacitance of the island. Let us denote the total ex- 
cess charge on the island by q — — ne. The electrostatic 
energy of the system can be written in a form where the 
charging energy is expressed as E c ^ = E c (n + n g ) 2 . Here 
E c = e 2 / (2Cs) and n g = V g C g /e is the gate charge. The 
change in energy for an electron tunneling to (+) or from 
(-) the island through junction i reads [7| 

I, (1) 



Eft* = ±2E c {n + n g ± 1/2) ± e(V 



where 



j- Ylj=i CjVj is the offset to the island po- 
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tential from the bias voltages. The <j> term can be inter- 
preted also as a shift to the gate position n g , and hence 
we neglect it in the following analysis where we consider 
the range of gate modulation of electric and heat current 
at fixed bias points. The standard expressions for the 
tunneling rate T and average heat flux Q as a function 
of the change in energy E% for each tunneling event are 
identical to those given in Ref. The charge number 
distribution p n is determined from the steady-state mas- 
ter equation In the analysis of the measured data, 
all quantities were evaluated numerically. The shape and 
position of the charge number distribution are crucial in 
describing the device properties in the Coulomb blockade 
regime fee In ^ E c , where Tn is the N island electron 
temperature. 

For an analytic investigation of Coulomb blockade ef- 
fects on heat transport near the optimal bias point where 
the maximum cooling power is achieved, we neglect the 
small probe current, and assume the two cooling junc- 
tions to have the same resistance Furthermore, we 



consider the extreme cases n g = (closed gate) and 
n g = —1/2 (open gate) in the regime E c /{k B T^) 3> 1, so 
that the charge number distribution p n is symmetric and 
narrow. Neglecting the population of the excited charge 
states, a symmetrically biased Coulomb blockaded SI- 
NIS with an open gate behaves as a regular SINIS whose 
electric and heat currents have been scaled by 1/2: At 
charge state n = (n — 1), an electron can enter (leave) 
the island through the junction biased at =fVds/2 with 
no increase in charging energy, while tunneling through 
the other junction is effectively blocked due to a charg- 
ing energy penalty of 2E C . We can therefore apply the 
known results valid in the absence of charging effects 
Q, namely that that the optimal cooler bias voltage is 
V£f ~ 2(A - 0.66fc B T N )/e, and the maximum cooling 
power for the heat transistor reads 
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With a closed gate, all tunneling events are affected by 
Coulomb blockade and the cooling power is minimized. 
The leading term of the cooling power corresponds to 
tunneling events occurring at the n — state, and em- 
ploying similar approximations as in Eq. l[2|). the cooling 
power accounting for the charging energy penalty of mag- 
nitude E c can be written 
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Comparison with numerical computations shows that the 
above approximative formulas are accurate to within 10% 
when E c /(k B T N ) > 2.5. 

In the experiment, we recorded the voltage over the 
current-biased probe junctions as the gate voltage was 
swept at a rate of ~ 1 Hz, while the voltage over 
the cooler junctions was swept significantly slower over 
a range of more than 2A/e. The sample parameters 
A = 230 /LtV and Cs = 1.5 fF were determined from a fit 
to the measured I-V characteristics and the gate mod- 
ulation amplitude of Vp ro be, respectively. We charted a 
number of bath temperatures ranging from 50 mK to 
560 mK. In Fig. [21 sweeps over cooler bias (Vbs) are 
shown at selected bath temperatures for both the refer- 
ence sample with vanishing charging energy [Fig. [Ha)] 
and for the heat transistor [Fig. [2jb)] . Due to very slow 
sweeping rate compared to the electron-phonon relax- 
ation time 0, each data point represents the system in 
thermal steady state. In the reference sample which be- 
haves as an ordinary SINIS cooler (i.e., with negligible 
charging energy), there is a one-to-one correspondence 
between the island electron temperature and the mea- 
sured probe voltage for a fixed current bias. The charac- 
teristics of the heat transistor are drastically different, as 
seen in Fig.[5Jb). In order to perform accurate thermom- 
etry on the heat transistor sample with varying V g and 
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FIG. 3: (color online) (a) Electron temperatures extracted 
from the data acquired at bath temperature 214 mK with the 
gate open (n g = 1/2, blue squares) and closed (n g = 0, red 
triangles), (b) The theoretical cooling power vs. — To 
for each data point of the left panel, and a linear fit to it. 
According to Eq. @, the slope is given by — EV. 



FIG. 2: (color online) Measured probe voltage as a function of 
the drain-source bias voltage for (a) the reference sample with 
negligibly small charging energy and (b) the heat transistor 
at different bath temperatures. For the heat transistor, the 
gate voltage was swept during the measurement, and hence 
the plot displays the full range of gate modulation at each bias 
point, whereas there was no observable gate modulation in the 
reference sample. Superimposed on panel (b) are theoretical 
curves for the two extremal gate positions. The solid blue 
curves are calculated assuming an electron-phonon heat load 
given by Eq. Q with the experimentally determined value 
£ = 2.3 x 10 9 WK -5 m~ 3 , and the dashed red curves given 
for reference show the expected behavior where the electron 
temperature is fixed at the bath tempeature (Tbath) over the 
full bias range. 



Vds 5 one needs to determine the temperature Tn by us- 
ing a complete model which takes into account the exact 
charge number distribution. 

Steady state is reached at a temperature where the 
cooling power is matched to the external heat load. The 
dominating source of heat load into the island electrons 
is electron-phonon heat flux (Qcp), which can be modeled 

by M 



Q cp = EV(T 5 - Tj|), 



(4) 



where E is a material parameter, V is the volume of N 
island, and To is the phonon temperature. Using bulk 
values from Ref. ll| for the Kapitza resistance between 
the island and substrate phonons, we estimate that the 
change in the temperature of the island phonons due to 
electron cooling is less than 10% of the observed drop in 
the electron temperature. The true Kapitza resistance 
for thin films is most likely even lower as one can as- 
sume a common phonon system for the film and sub- 



strate |12ll. and thus we neglect lattice cooling unlike in 
Ref. [13j, and use To = 7t> a th- In our modeling, we fur- 
ther ignore heating of the superconducting reservoirs due 
to hot quasiparticles extracted from the island, whose ef- 
fect is supposedly weak due to relatively high junction 
resistances and moderate cooling power. Consequently, 
we take Ts = Tbath- On the other hand, at sub-100 mK 
bath temperatures, electrons are overheated due to noise 
via the electrical leads, and accurate comparison with the 
theoretical model becomes difficult. 

In analysis of the measurement data, we concentrate 
on characterizing the cooling properties of the device at 
the extremal gate positions. These have to be identi- 
fied as the local minima and maxima of either V pro hc or 
the drain-source current as the gate voltage is swept. In 
Fig. El^a), we present the electron temperatures corre- 
sponding to open and closed gate positions derived from 
the experimental data at bath temperature 214 mK. Sim- 
ilar results were obtained at other bath temperatures 
as well, but only in the range 170 - 250 mK we have 
both accurate thermometry and significant temperature 
reduction. The temperature Tn was evaluated numeri- 
cally for each data point using the experimentally deter- 
mined sample parameters, the observed values of lp ro be 
and Vds, and the set value of probe current, which was 
3.2 pA in this measurement. The electron temperatures 
obtained from the model for small bias voltages Vus, 
for which no cooling is expected, differ less than 10 mK 
from the bath temperature for both open and closed gate, 
demonstrating a good agreement between the experiment 
and theory, and supporting the reliability of this method 
of thermometry. With open gate, the transistor should 
function almost like a regular SINIS cooler, and cooling 
peaks are indeed found at bias voltages Vds — ±2A/e. 
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Increasing the bias voltage over ~ 2A/e brings the device 
rapidly into the ohmic heating regime, which is observed 
as a sharp rise of the electron temperature. For the closed 
gate configuration, the cooling peaks are weaker and oc- 
cur at a slightly larger bias voltage, as tunneling electrons 
have to overcome both the Coulomb and the supercon- 
ducting gap. 

A quantitative analysis of the observed temperature 
drops is presented in Fig. [3^b). The plot shows the nu- 
merically calculated total cooling power versus T^ — Tq , 
using the Tn values from Fig.[3Ja). According to the heat 
balance equation Q = Q cp , the dependence should be lin- 
ear with the slope given by — EV. Agreement with the 
model is excellent throughout the whole data range from 
cooling (Q > 0) to resistive heating (Q < 0). Through 
a linear fit we obtain £ = (2.3 ± 0.3) x 10 9 WK- 5 nr 3 , 
the main source of uncertainty being the effective island 
volume. The literature value E = 2 x 10 9 WK -5 m~ 3 
[3] is within the experimental error. Estimates for E 
obtained from data at other bath temperatures are con- 
sistent with the above result, demonstrating the stability 
of the method. Furthermore, the observed on/off ratio 
of cooling power upon gate variation at the optimal bias 
point was as large as 3.2. 

Figure 12(b) displays the V^robe curves (blue lines) cal- 
culated using the model which accounts for changes in 
the island electron temperature with the experimentally 
obtained value E = 2.3 x 10 9 WK _5 m~ 3 . For reference, 
we also show curves with the electron temperature fixed 
to Tbath (red dashed lines), demonstrating which of the 
observed features should be attributed to changes in tem- 
perature. All the qualitative features of Vp ro be are repro- 
duced, including the gentle bumps appearing at higher 
temperatures, which are not indicative of cooling, but 
of changes in the width of the charge number distribu- 
tion. These features can be observed in V pYO h e when 
E c < &bTn, Tn ~ Ts and when Vp TO be is small. As ex- 
pected, significant differences between the cooled Tn and 
Tn = Tbath curves appear only at low temperatures and 
with open gate, whereas with closed gate, features due to 
change in the island electron temperature remain always 
small. The measured amplitude of gate modulation at 
Tbath = 137 mK is slightly smaller than expected. The 
difference can be attributed to heat load from electrical 
noise, to which the measurement is increasingly sensitive 
at low temperatures. The smallness of the cooling power 
of the probe current can be verified by the almost vanish- 
ing difference between the theoretical curves for constant 
and variable Tn near Vus = 0. At high bath temper- 
atures, the measured amplitude of the gate modulation 
and the shape of the envelope curves are as expected, 
although there are unexplained offsets of the order of 
10 /j,V in the measured data at different bath tempera- 



tures as compared to theoretical curves. We have also 
verified that the theoretical model reproduces accurately 
the measured cooler junction I-V curves at all bath tem- 
peratures (data not shown), which further validates the 
accuracy of the model and the fitted sample parameters. 

In conclusion, we have investigated experimentally and 
theoretically the thermal DC properties of a Coulomb- 
blockaded SINIS structure. We have demonstrated that 
the system operates as a mesoscopic heat transistor, i. e., 
it shows a gate modulation of heat flux. In particular, 
cooling power modulation as large as 70% of its max- 
imum has been observed. The agreement between the 
experiment and a theory based on Coulomb-blockaded 
single-electron tunneling is compelling. Furthermore, we 
have also proved that an NIS probe can be used as a 
sensitive thermometer even in the presence of charging 
effects, provided that the charge number distribution on 
the normal metal is carefully modeled. The influence of 
this distribution can be reduced up to large extent by se- 
lecting the bias point properly. Finally, our work might 
provide guidelines for the investigation of heat transport 
at the mesoscopic scale, e.g., to either avoid or make use 
of the effects arising from single-electron charging in elec- 
tronic refrigeration. 
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